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We present results on the mass of the nucleon and the A using two dynamical degenerate twisted 
mass quarks. The evaluation is performed at four quark masses corresponding to a pion mass 
in the range of 690-300 MeV on lattices of size 2.1 fm and 2.7 fm. We check for cutoff effects 
by evaluating these baryon masses on lattices of spatial size 2.1 fm with lattice spacings fl(j3 = 
3.9) = 0.0855(6) fm and fl(j3 = 4.05) = 0.0666(6) fm, determined from the pion sector and 
find them to be within our statistical errors. Lattice results are extrapolated to the physical limit 
using continuum chiral perturbation theory. The nucleon mass at the physical point provides a 
determination of the lattice spacing. Using heavy baryon chiral perturbation theory at 6'[p^) we 
find a(j3 = 3.9) = 0.0879(12) fm, with a systematic error due to the chiral extrapolation estimated 
to be about the same as the statistical error This value of the lattice spacing is in good agreement 
with the value determined from the pion sector. We check for isospin breaking in the A-system. 
We find that A++-^ and A+ '' are almost degenerate pointing to small flavor violating effects. 
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1. Introduction 

Twisted mass fermions provide a promising formulation of lattice QCD that allows for auto- 
matic i^{a) improvement, infrared regularization of small eigenvalues and fast dynamical simula- 
tions [|I|]. We use the tree-level Symanzik improved gauge action and work at maximal twist to 
realize ^(a) -improvement. Recent results obtained in the pion sector give an accurate evaluation 
of the low energy constants Tj, and /4 ||], which lead to the most accurate determination of the 
S-wave Tin scattering lengths [Q]. In this work we study the light baryon sector. 

The fermionic action for two degenerate flavors of quarks in twisted mass QCD is given by 

Sf = a'^Y.'^(^){^w[U]+mo + i^Y5Z^)\j/{x) (1.1) 

X 

with Dw [U ] the massless Dirac operator, niQ the bare untwisted quark mass and /i the bare twisted 



mass. The twisted mass term in the fermion action of Eq. (1.1) breaks isospin symmetry since the u 



and d-quarks differ by having opposite signs for the /i-term. This isospin breaking is a cutoff effect 
of G{cP-\ However the up- and down-propagators satisfy Guix^y) = 75G^(3',x)75, which means that 
two-point correlators are equal with their hermitian conjugate with u- and d-quarks interchanged. 
Since the masses are computed from real correlators this leads to the following pairs being de- 
generate: 71+ and 71^, proton and neutron and A++(A+) and A^(A"). A theoretical analysis 
shows that potentially large &{a^) effects that appear in the 7r"-mass are suppressed in all other 
quantities. Calculation of the mass of tt", which requires the evaluation of disconnected diagrams, 
has been carried out confirming large ^(a^) -effects. In the baryon sector we can study isospin 
breaking by evaluating the mass difference between A++(A^) and A+(A°). Since no disconnected 
contributions enter we can obtain an accurate evaluation of isospin splitting and its dependence on 
the lattice spacing. 

2. Lattice techniques 

The parameters of the calculation are collected in Table [T[ They span a pion mass range from 
300-690 MeV. At a pion mass of about 300 MeV we have simulations for lattices of spatial size, 
Ls = 2.1 fm and Lg = 2.7 fm at j8 = 3.9 allowing to check finite size effects. We provide a first 
check of finite a-effects by comparing results at /3 = 3.9 and j3 = 4.05. 

The masses of the nucleon and the A's are extracted from two-point correlators using the 
standard interpolating fields, which for the proton, the A++ and A+, are given by 



Jp = EabciulCysdhjuc, Ja++ = eabc{ulC'Y^Ub)uc 



1 

•^A+ — ^abc 



2[u^Cy^db)uc+ {u^CY^Ub)dc 



(2.1) 



Local interpolating fields are not optimal for suppressing excited state contributions. We instead 
apply Gaussian smearing to each quark field, q{x,t): q^™'^'^^{x,t) = Y.yF{x,y;U{t))q{y,t) using the 
gauge invariant smearing function 

Fix,y;U{t)) = {\ + aH)"{x,r,U{t)), (2.2) 
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= 3.9, a 


= 0.0855(6) imimmfn % 




X 48, L, 


= 2.1 fm 






0.0040 0.0064 0.0085 


0.010 








(GeV) 


0.3131(16) 0.3903(9) 0.4470(12) 


0.4839(12) 


32^ X 64, Ls 


= 2.7 fm 






0.004 










(GeV) 


0.3082(55) 










= 4.05, a 


= 0.0666(6) fmfrom/;, ^ 




32^ X 64, L, 


= 2.1 fm 




At 


0.0030 0.0060 0.0080 


0.010 








(GeV) 


0.3070(18) 0.4236(18) 0.4884(15) 


0.6881(18) 



Table 1: The parameters of our calculation. 



constructed from the hopping matrix, H{x,y;U{t)) = LLi (^^K'^'O^x.y-f + ?//(x — /,?)5x,y+,J . 

The parameters a and n are varied so that the root mean square (r.m.s) radius obtained using the 
proton interpolating field is in the range of 0.3-0.4 fm. In Fig. |l| we show lines of constant r.m.s 
radius as we vary a and n. The larger the n the more time consuming is the smearing procedure. 
On the other hand, for a ~ 1, increasing further a does not reduce n significantly. Therefore, 
we choose a value of a large enough so that the weak ^-dependence sets in, and we adjust n to 
obtain the required value of the r.m.s radius. We consider two sets for these parameters giving r.m.s 
radius 0.31 fm and 0.39 fm, as shown in Fig. |l]. In Fig. ^, we show the nucleon effective mass, 

= — log(C(f)/C(f — 1)) with C{t) the nucleon correlator, for 10 configurations at jS = 3.9 and 
/i = 0.0085. For the optimization of the parameters we apply Gaussian smearing at the source, 
whereas for the sink we use local interpolating fields so that no additional inversions are needed 
when we change a and n. As can be seen, for both sets of smearing parameters, the excited state 
contributions are suppressed with the set a = 4, « = 50 producing a plateau a couple of time slices 
earlier. If, in addition, we apply APE smearing to the spatial links that enter the hopping matrix, 
then gauge noise is reduced resulting in a better identification of the plateau. Therefore for all 
computations at j8 = 3.9 we use Gaussian smearing with a = 4 and n = 50. We apply smearing at 
the source and compute the mass using both local (LS) and smeared sink (SS). For j8 = 4.05 we 
readjust the parameters so that the nucleon r.m.s radius is still about 0.39 fm, obtaining a = 4 and 
n = 70. In all cases we apply APE smearing to the gauge links that are used in F(x,y;U [t)). We 
note that Gaussian smearing is very effective as compared to, for example, fuzzing on links joining 
quarks at different sites. 

The nucleon effective masses obtained using correlators with smeared source and local or 
smeared sink for the four /i-values at j3 = 3.9 are shown in Fig. ^, where we average over the 
proton and neutron correlators. In Fig. ^ we show, for the same /x-values, the A effective masses 
after averaging the correlators obtained using smeared source and sink over the degenerate pairs 
A++, A^ and A+, A° . As can be seen, the quality of the plateaus in the nucleon case is better than 
in case of the A. This explains why results on the A mass have larger errors requiring more statistics 
for a reliable determination. The errors are evaluated using jackknife and the F-method ^ to check 
consistency. In all the figures we show the errors obtained with the latter method. 

3. Results 

We show our results for the nucleon mass as a function of in Fig. |5|, where we use a = 
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Figure 1: Lines of constant r.m.s radius as 
function of the smearing parameters a and 
n. The asterisk shows the values a — 2.9, 
n = 3Q and the cross a — 4.0, n — 50. 
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Figure 2: versus time separation both in lattice units. 
Crosses show results using local sink and source (LL), cir- 
cles (asterisks) using Gaussian smearing at the sink (SL) with 
a — 2.9 and n = 30 ( a = 4 and n ~ 50), and filled triangles 
with a = 4 and n = 50 and APE smearing. The dashed line 
is the plateau value when APE smearing is used. 
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Figure 3: Nucleon effective mass (LS: asterisks, SS: Figure 4: A++'^ (asterisks) and A^'" (open triangles) 

open triangles) for j3 = 3.9 versus time separation in effective masses for fi —3.9 versus time separation in 

lattice units. lattice units. 



0.0855 at j8 = 3.9 and a = 0.0666 at j3 = 4.05, determined from [|], to convert lattice results 
to physical units. As can be seen, the results at these two j3 -values show good scaling pointing 
to small cutoff effects. For the three larger pion masses nijiLs > 4, whereas for the smallest value 
m;r(/i = 0.004)Li ~ 3.2. Applying the resummed Luscher formula to the nucleon mass and using 
the knowledge of the nN scattering amplitude to iff{p^) and ^{p^) it was shown that, for ~ 2 fm 
and nijc ~ 300 MeV, the volume corrections are small being estimated to be about (3-5)% |^]. 
We calculate the nucleon mass increasing the spatial length of the lattice from 2.1 fm to 2.7 fm 
so that mjc{lJ. = 0.004)Li ~ 4.3. If Ahtat = mi^{Ls/a = 24) —mi^{Ls/a = 32) then we find that 
l^mf^ jmt^iLsl a = 32) = 0.01 ±0.02 at our smallest quark mass, i.e. consistent with zero within 
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our statistical error but also within the estimated error range of Ref. [^]. In Fig. ^ we include, for 
comparison, results obtained with dynamical staggered fermions from Ref. As can be seen, the 
results using these two formulations are consistent with each other. 

In Fig. ^ we show our results for the mass difference between the averaged mass of the pairs 
A++, and A+, A". As can be seen, the splitting is consistent with zero, indicating that isospin 
breaking in the A system is small. 
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Figure 5: The nucleon mass as a function of for 
j3 = 3.9 on a lattice of size 24^ x 48 (filled triangles) 
and on a lattice of size 32^ x 64 (open triangles). Re- 
sults at j3 = 4.05 are shown with the stars. The phys- 
ical nucleon mass is shown with the asterisk. Results 
with dynamical staggered fermions for A^/r = 2 + 1 
(filled circles) and Nf —2 (open circle) on a lattice 
of size 20^ X 64 with a = 0.125 fm are from Ref. [Ml. 



Figure 6: The mass splitting between A+'*' and A++-^ 
as a function of m\ both in lattice units. The notation 
is the same as in Fig. 5. 



Having checked that \mi^ at the smallest pion mass is consistent with zero within our statistical 
errors and that cut-off effects are small, we use, in what follows, continuum chiral perturbation 
theory in an infinite volume to perform the chiral extrapolation to the physical point. The leading 
one-loop result in heavy baryon chiral perturbation theory (HB;^PT) |^ is well known: 



rriN 



■■ nil. 



■4ciml- 



3gl 



■2 '"^ 



(3.1) 



with m^, the nucleon mass at the chiral limit, and ci treated as fit parameters. We find that this 
) result provides a very good fit to our lattice data at j3 = 3.9, yielding m% = 0.875(10) GeV 
and c\ = —1.23(2) GeV^^ with ;f/d.o.f. = 0.2. In this determination we use a = 0.0855 and 
results obtained on both lattice volumes. The value extracted for c\ can be compared to the value 
c\ = —0.9 ±0.5 GeV^^ extracted from various partial wave analyses of elastic TiN scattering data 
for the TiA/^-sigma term. We would like to stress that, despite the fact that the physical point is not 
included in the fit as customary done in other chiral extrapolations of lattice data, the nucleon mass 
that we find at the physical pion mass is 0.955(10) GeV. Given that the error is only statistical, the 
fact that this value is so close to the experimental value is very satisfactory. Chiral corrections to 
the nucleon mass are known to ^{p^) within several expansion schemes. In this work to ^{p^) we 
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use the results obtained in HB;^PT with dimensional regularization 1 10] and in the so called small 



scale expansion (SSE) [|11[]. HB;^PT with dimensional regularization is in agreement with covariant 
baryon ;^PT with infrared regularization up to a recoil term that is of no numerical significance. In 
SSE the A-degrees of freedom are explicitly included in covariant baryon ;^PT by introducing as 
an additional counting parameter the A-nucleon mass splitting, A = mA — rriN, taking 0'{A/mN) ~ 
0'{mjt/mN). A different counting scheme, known as 5-scheme, takes A/rriN ~ ^{S) and mji/niN ~ 
ff{5^) [|I2|]. Using the 5— scheme in a covariant chiral expansion to order p^, p^/A one obtains 
an expansion that has a similar form for the nucleon and A mass. Here we use the variant of the 
5-scheme that includes the TiA-loop and adds the fourth order term C2m^ as an estimate of higher 
order effects, since the complete fourth order result is not available. The parameter C2 is to be 
determined from the lattice data. The fits using these different formulations are shown in Fig. ^. 
All formulations provide a good description of the lattice results and yield a nucleon mass at the 
physical point that is close to the experimental value. The physical nucleon mass is not including 
in the fits. We can use these chiral expansions to fix the lattice spacing using the nucleon mass at 
the physical point and compare with the value determined from the pion sector. The results of the 
fits in HB;^PT to ff{p^) and ^(/?^) are shown in Fig. ^. Using the leading one-loop result we find 
a = 0.0879(12) fm, whereas to ^{p'^) we obtain a = 0.0883(9) fm. Both SSE and the 5-scheme, 
which include explicitly A-degrees of freedom, yield values that are consistent with those obtained 
in HB;^PT. The variation in the value of a in the different chiral extrapolation schemes gives an 
estimation of the systematic error involved in the chiral extrapolation. A proper determination of 
the systematic error is in progress. 
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Figure 7: Chiral fits to the nucleon mass using 
a = 0.0855. The physical point shown by the as- 
terisk is not included in the fits. 



Figure 8: Chiral fits using HBjPT to ^(p^) and 
0{p'^) determining a from the nucleon mass. 



The leading one-loop HB;^PT result in the case of the A mass has the same form as that for 
the nucleon mass given in Eq. ( ^ ) with —>■ and ci —>■ cia- Assuming SU(6) symmetry, 
the one-loop contribution has the same numerical value as in the nucleon case. It is useful to 
chirally extrapolate the A mass to see how close current results are to A(1232) taking a = 0.0879 
as determined from the nucleon mass. In Fig. we show the resulting fit. 
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Lattice data on the A++'^ mass, chirally ex- 
trapolated using tlie ^(7?^)-result in HB;^PT, 
yield, at the physical point, m^++,- = 
1.265(26) GeV, consistent with the resonant 
A mass. A similar chiral fit to the A+''' mass 
yields a curve that lies above the physical 
point but with an overall statistical error band 
that overlaps the one obtained from the chiral 

"m .05 .10 .15 .20 .25 .30 fit to the A++-" mass. 

(GeV^) 

Figure 9: Chiral fits to the A++'^ and mass using HBjFT to ffip^) with a set from the nucleon mass. 

4. Conclusions 

We have shown that twisted mass QCD yields accurate results on the nucleon mass close to 
the chiral regime. The quality of the results for pion masses in the range of 300-500 MeV allows a 
chiral extrapolation using heavy baryon chiral perturbation theory to ff{p^). The nucleon mass at 
the physical point provides a good physical quantity for setting the scale. Using the leading one- 
loop result in HB;^PT we find a{p = 3.9) = 0.0879(12) fm. Comparing this value to the results 
obtained using higher order terms in the chiral expansion, gives a first estimate of the systematic 
uncertainty due to the chiral extrapolation, that is of the same order of magnitude as the statistical 
error. Within this estimated uncertainty of the chiral extrapolation, the value we find for a(j8 = 3.9) 
at leading order in HB;^PT is consistent with the value determined from fji. The mass splitting in 
the A isospin multiplets calculated with two lattice spacings on two volumes is consistent with zero 
showing that isospin breaking effects are not severe in this channel. 
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